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for (y=0; y<M+3; ++y) { 
for tx=0; x<N+5; ++x) { 

if ({x-3)>=l && (x-5)<=N-2 {y-2)>=l (y-3)<^M-2) { 
if t(x-5)>=l && (y-3)>=l) I 
if (out_coTnpute =5= 255) { 

if (comp_edge_pixels [ ( (x-4) %3} *3 + (y-2) %3] <comp_edge_middle) out_compute=0; 

if (comp_edge_jpixels [ ( (x-4) %3) *3 + (y-4) %3] <comp_edge_Tniddle) out_compute=0; 

if {comp_edge_pixels [ ( (x-5)%3) *3+<y-4) %3] <comp_edge_middle) out_compute=0 ; 

if ({x-3)<=N-2 && (y-2}<=(M-2)) { 
maxd i f f _coTnpu t e = 

maxl3 (abs (gauss_xy_pixels [ ( (x-2) %3) *3+ (y-1) %3] 

- gauss_xy_tniddle) , maxdif f_compute) ; 

maxd i f f _compu t e = 

inaxl3 (abs (gauss_xy_pixels [ { (x-2) %3) *3+ (y-3) %3l 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute = 

inaxl3 (abs (gauss_xy_pixels [ ( (x-3 ) %3) *3+ (y-3) %3] 

- gauss_xy_middle) , maxdif f_compute) ; 
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for (y=0; y<M+3; ++y) { 
for (x=0; x<N+5; +-t-x) { 

if (x>=4 && x<=N+3 y>=3 y<=M+l) { 
if{{x-5)>-l && (y-3)>*l} { 
if (out_compute 255) { 

csexmin4mod3x3 = ({x-4)%3)*3; 
cseymin4mod3 = (y-4)%3; 

if (conip_edge_pixels [csexmin4Tnod3x3 + (y-2 } %3] <comp_edge_middle) out_cotnpute=0 ; 

i f {comp_edge_j>ixels [csextnin4mod3x3 +cseyTnin4mod3 3 <con^_edgejmiddle) out_compute=0 ; 

if tcomp_edge_pixels [ ( {x-5) %3 ) *3 +cseymin4mod33 <coinp_edge_middle) out_compute=0; 

if ((x-3)<=N-2 && (y-2)<=(M-2) ) { 
csexmin2mod3x3 ^ {(x-2)%3}*3; 
cseymodS = y%3; /* = (y-3)%3 */ 
maxdif f_compute = 

maxl3 (abs (gauss_xyjpixels [csexmin2inod3x3 + (y-1) %3] 

- gauss_xyjnEiiddle) , maxdif f_coTnpute) ; 

maxdif f_cornpute * 

maxl3 (abs (gauss_xyjpixels [csexmin2mod3x3 + cseymod3] 

- gauss_xy_iniddle) , maxdif f_conpute) ; 

maxdif f__corT^ute 

maxl3 (abs (gaus_xy_pixels { (x%3) *3 + cseytnod3] 

- gauss_xy_middle) , maxdif f_compute) ; 
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for (y=0; y<M+3; ++y) { 
cseymod3 = y%3; 
cseyminltaod3 (y-l)%3; 

cseymin2mod3 = (y-2)%3; ^ 
cseymin4mod3 s= (y-4)%3; 
for (x-0; x<N+5; ++x) { 

if {x>=4 && x<=N+3 && y>=3 && y<=M+i) { 
if((x-5)>=l && (y-3)>=l) \ 
if (out_compute ==255) { 

csexmin4mo<a3x3 = {{x-4)%3)*3; 

if {coinp_edge_pixels Ecsextnin4mod3x3 +cseymin2mod3l <comp_edge_middle) out_coTrvpute=0; 

if (coTnp_edge_pixels Ecsexmin4Tnod3x3 +cseymin4mod3] <comp_edge_middle) out_compute=0; 
if (coinp_edge_pixels t ( (x-5) %3 ) *3 +cseymin4mod3l <comp_edge_iniddle) out_compute=0 ; 

if ((x-3)<=N-2 && (y-2)<=(M-2) ) { 
csexmin2mod3x3 = ((x-2)%3)*3; 
maxdiff^compute = 

maxl3 {abs (gauss_xy_pixels Ecsexmin2mod3x3 + cseyminlTnod3] 

- gauss_xy_iniddle) , maxdif f_compute) ; 

maxdiff_compute = 

maxl3 (abs (gauss_xy_jpixels [csexmin2inod3x3 + cseymod3} 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute = 

maxl3 (abs (gaus_xy_j>ixels t (x%3) *3 + cseymod3] 

- gauss_xy__iniddle) , maxdif f_coinpute) ; 
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cseymod3 = -1; 

for (y=0; y<M+3; ++y) { 
cseyminlmod3 = cseyinod3; 
cseymod3 = y%3; 
cseymin2mod3 = 3-cseyinod3-cseyininlmod3 ; 
for (x«0; x<N+5; ++x) { 

if (x>=4 && x<=N+3 && y>=3 && y<=M+l) { 
ift(x-5)>=^l && (y-3)>=l) { 
if (out_compute == 25S) { 

csexmin4mod3x3 = ( {x- 4 ) %3 ) *3 ; 

if ( coinp_edge_pixels [csexmin4mod3x3 +cseymin2mod3] <coTnp_edge_middle) out_compute=0 ; 

if (comp_edge_pixels [csexmin4mod3x3 +cseyminlmod3] <comp_edge_middle) out_coTnpute=0 ; 
if (coTnp_edge_j)ixels [ ( (x-5) %3) *3 +cseyminlmod3l <comp_edge_middle) out_compute=0 ; 



if ({x-3)<=N-2 && {y-2)<=(M-2)) { 
csexmin2Tnod3x3 = ( (x-2) %3) *3 ; 
maxdif f_compute = 

iTiaxi3 (abs (gauss_xy_jpixels [csexinin2mod3x3 + cseyminlmod33 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_conipiite - 

maxl3 (abs {gauss^xyjpixels [csexmin2mod3x3 + cseymodS] 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_conipute « 

maxl3 tabs (gaus_xyjpixels E (x%3) *3 + cseymod33 

- gauss_xy_middle) , maxdif f^compute) ; 
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cseymodS = -1; 

for {y=0; y<M+3; ++y) { 

cseyminlmodB = cseymodS; 

cseymodS = y%3; 

cseymin2raod3 = 3-cseymod3-cseyminlmod3; 
for (x=0; x<N+5; ++x) { 

if (xi=4 && x<=N+3 £c& y>=3 && y<=M+l) { 
csexmod3x3 = {x%3)*3; 
csexmin2Tnod3x3 = ((x-2)%3)*3; 
csexmin4inod3x3 = ((x-4)%3)*3; 
csexmin5niod3x3 = ((x-5)%3)*3; 
if({x-5)>=l &&r (y-3)>=l) f 
if (out_compute == 255) { 

if {coTnp_edge_pixels tcsexmiii4mod3x3 +cseymin2mod3]<comp_edge_niiddle) out_compute=0; 

if {coTnp_edge_pixels [csexmin4mod3x3 +cseyT[iinlniod3] <coinp_edge_middle) out_compute=0; 

if (coinp_edge_pixels tcsexmin5mod3x3 +cseyminlmod3] <comp_edge_middle) out_compute=0; 

if ((x-3)<=N-2 && (y-2)<=(M-2)) { 
maxdif f_corapute = 

maxl3 (abs (gauss_xy_pixels [csexTtiin2niod3x3 + cseyminlmod3] 

- gauss_xy_iiiiddle) , maxdif f_compi2te) ; 

maxdif f_compute = 

maxl3 (abs (gauss_xy_pixels [csexmin2mod3x3 + cseymod3] 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_coinpute = 

maxl3 (abs {gaus_xy_j)ixels [csexmod3x3 + cseymod33 

- gauss_xy_middle) , maxdif f_compute) ; 
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cseymod3 = -1; 

for (y=0; y<M+3; ++y) { 

cseyminlmod3 = cseymodB; 

cseymodS = y%3; 

cseyinin2inod3 = 3-cseymod3 -cseyTninlmod3 ; 
for (x=0; x<N+5; ++x) { 

csexinod3x3 = {x%3)*3; 

csexmin2mod3x3 = ((x-2)%3)*3; 

csexmin4mod3x3 = ((x-4)%3)*3; 

csexmin5Tnod3x3 = ((x-5)%3)*3; 

if (x>=4 && x<=N+3 && y>=3 && y<=M+l) { 
if{(x-5)>=l && (y-3)>=l) { 
if {out_contpute == 255) { 

if {coinp_edge_pixels tcsexmin4mod3x3 +cseymin2mod3] <comp_edge_middle) out_compute=0 ; 

if (comp_edge_pixels tcsexmin4mod3x3 +cseyminlmod3} <comp_edge_middle) out_compute=0; 

if {comp_edge_pixels [csexmin5mod3x3 +cseyminlmod3} <comp_edge_middle) out_compute=0 ; 



if ((x-3)<=N-2 && (y-2)<=={M-2)) { 
maxdif f_coTnpute = 

maxl3 (abs (gauss_xy_pixels [csextttin2niod3x3 + cseyminlniodS] 

- gauss_xy_middle) , maxdif f^compute) ; 

maxdif f__compute = 

maxl3 (abs (gauss_xy_pixels [csexmin2mod3x3 + cseymod3] 

- gauss_xy_middle} , maxdif f_coTnpute) ; 

maxdif f_cornpute = 

niaxX3(abs {gaus_xy_j>ixels [csexinod3x3 + cseymod3j 

- gauss_xy_middle) , maxdif f_compute) ; 
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cseymod3 = - 1 ; 

for {y=0; y<M+3; 4+y) { 

cseyrainlmodB = cseymod3; 

cseymod3 = y%3; 

cseymin2niod3 = 3-cseymod3-cseyminlmod3 ; 

csexmod3x3=-3; 

for {x=0; x<N+5; ++x) { 

csexiiiinlmod3x3 = csexmod3x3 ; 

csexmod3x3 = (x%3)*3; 

csexmin2inod3x3 = 9-csexmod3x3-csexminlmod3x3; 

if tx>=4 && x<==N+3 && y>=3 && y<=M+l) { 
if ( (x-5) >=1 && [y-3)>=l) { 
if (otit_coTnpute == 255) { 

if (comp_edgejpixels tcsexminlmod3x3 +cseymin2mod3] <coTnp_edge_middle) out_coTnpute=0 

if (comp_edgejpixels [csexniinlmod3x3 +cseyininltnod33 <comp_edge_middle) out_coTnpute=0 

if (comp_edge_pixels [csexmin2mod3x3 +cseyminlmod3] <comp_edge_middle) out_coinpute=0 



if ((x-3)<=N-2 && {y-2)<={M-2)) { 
maxdif f_compute = 

maxl3 tabs {gauss_xy_pixels [csexmin.2inod3x3 + cseyminlmod3] 

- gauss_xy_Tniddle) , maxdif f_compute) ; 

maxdif f_compute = 

Tnaxl3 (abs (gauss_xy_pixels tcsexmin2mod3x3 + cseymod3] 

- gauss_xy_middle) , maxdif f_compute) ; 

maxdif f_compute = 

maxl3 (abs (gaus_xy_pixels [csexmod3x3 + cseymod3] 

- gau s s_xy_mi dd 1 e ) , maxd i f f _comp\i t e ) ; 
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cseymod3 = -1; 

for (y=0; y<M+3; ++y) { 

cseyminlTnod3 = cseymodS; 

cseyniod3 = y%3 ; 

cseymin2mod3 = 3-cseymod3-cseyminlmod3 ; 
csexmod3x3=-3; 
cseyminlTnod2 = (y-l)%2; 
cseymod2=l-cseyminlmod2 ; 
for (x=0; X<N+5; ++x) { 

csexminlmod3x3 = csexmod3x3; 

csexmod3x3 - (x%3)*3; 

csexmin2mod3x3 = 9-csexmod3x3-csexminlTnod3x3 ; 
csexminlx2 = (x-l)*2; 
csexmin3x2 = csexminlx2-4 ; 

if (x>=3 && x<N+3 y>=2 y<M+2) 
tmparray [ (csexmin3x2+ cseymod2) %160 

+ (csexmin3x2+ cseymod2) / 160*256 + 96} 
= comp__edgej)ixels Icsexmod3x3 

+ cseyTnin2inod3] = maxdif f _compute 

if (x>= 1 x<N+l && y>=l y<=M) 

tmparray [ (csexminlx2 + cseyminltnod2) %64 

+ (csexminlx2 + cseyminlmod2) /64*2563 
- gauss_xy__pixels Ecsexminlmod3x3 

-t- cseyminlmod3] = gauss_xy_coTnput 



} 




cseyTnod3=-l; 

for (y=0; y<M+3; ++y) { 

cseyminlmod3=cseymod3 ; 

cseymod3 ++ ; 

if(cseymod3 >= 3) {cseymod3 -= 3;} 

cseymin2mod3 = 3-cseymod3-cseyminlmod3 ; 

cseytninlmod2 = (y-l)&l; 

cseyTnod2=5^1-cseyininlmod2 ; 

csexmod3x3=i -3; 

csexx2modl 6 0_l_2=cseymod2 - 8 ; 

csexx2divl60_l_2=0; 

for (x=0; x<N+5; ++x) { 

csexminlmod3x3=csexmod3x3 ; 

csexmod3 ++ ; 

if(csexmod3 >= 3){csexmod3 -= 3;} 
csexmod3x3=csexTnod3*3 ; 

csexmin2mod3x3 = 9-csexmod3x3-csexminlmod3x3; 
csexx2inodl60_l_2+=2 ; 
csexminlx2 = {x-l)-*2; 
csexmin3x2 = csexminlx2-4 ; 

if {csexx2inodl60_l_2>=160) {csexx2modl60_l_2-=160;csexx2divl6 0_l_2++; } 

if (x>=3 && x<N+3 y>=2 && y<M+2) 

tmparray tcsexx2inodl60_l_2 + csexx2divl60_l_2*256 + 96] 
= comp__edge_pixels [csexmod3x3 

+ cseymin2mod3] = maxdif f_coTnpute; 

if (x>== 1 x<N+l && y>=l && y<=M) 

tmparray t { (csexminlx2 + cseyminlmod2) &63 ) 

+ ( (csexminlx2 + cseyminlmod2) >>6) *256] 
= gauss_xy_pixels [csexminlmod3x3 

+ cseyminlmod3l = gauss_xy_compute; 
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Consider expressions in 2 different 
conditionals . 
(1) The 2 conditionals ranges are over lagging ? 




. . The expressions 

shows conanon factors ? 




No CH for these expressions 




(4)rangel + range 2 > whole range ? 




No CH for these expressions 



(5) Evaluate: 

- the ranges compared to the whole range (kl and k2) 

- the cost of expressions in each branch (cl and c2) 

- and the cost after a potential CH {C_new) 

- degrees of similarity {S « (cl+c2-C_new) /C_new ) 

C1+C2 < {1+S)*{cl*kl + C2*k2) ? 



Do CH across 
conditional for 
these expressions 



No CH for these 
ej^ressions 





Do CH across 
conditional for 
these expressions 

S=0 => (l>/no 



No CH for these 
expressions 
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